Improved measurement of time-dependent CP violation in B J/ip-K decays 
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We report improved measurements of time-dependent CP violation parameters for _B''(B°) —f 
J/ip-K° decay. This analysis is based on 535 million BB pairs accumulated at the T(4S') resonance 
with the Belle detector at the KEKB asymmetric-energy e"^e~ collider. From the distribution of 
proper time intervals between the two B decays, we obtain the following CP violation parameters 

Sj/i,.^o = -0.65 ± 0.21(stat) ± 0.05(syst) and 
Az/v-ttO = +0.08±0.16(stat) ±0.05(syst), 

which are consistent with Standard Model expectations. 
PACS numbers; 12.15.Hh, 13.25.Hw 



The Kobayashi-Maskawa (KM) quark-mixing ma- 
trix [ij has an irreducible complex phase that gives rise 
to CP- violating asymmetries in the time-dependent rates 
of and decays into a common CP eigenstate such 
as J/V'7r° 0. In the decay chain T(4S)^ B^W ^ 
( J/^/; 7r°)/tag, where one of the B mesons decays at time 
tcp to the final state J/ipn^ and the other decays at 
time itag to a final state /tag that distinguishes between 
B° and B^ , the decay rate has a time dependence given 

by i 



p-|At|/r^o 
V{At) = —^ ^1 + q- 



sin(ATOdAi) 



where Tgo is the neutral B lifetime, Arud is the mass 
difference between the two neutral B mass eigenstates, 
At = tcp — ttag, and the 6- flavor charge q — +1 (—1) 
when the tagging B meson is a B^{B'^). The CP viola- 
tion parameters Sj /^t^o and Aj/^irO are given by 
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where A is a complex parameter that depends on both 
the B'^B'^ mixing and the amplitudes for 5° and de- 
cay to J/ipTT^. In the Standard Model (SM), |A| is, to a 
good approximation, equal to the absolute value of the 
ratio of the J/il) tt^ to B^ J /ip tt" decay am- 

plitudes. At the quark level, the B^ J/i/jtv'^ decay 
proceeds via a 6 — > ccd transition. In this decay, the 



tree amplitude is CKM-suppressed. Since the tree ampli- 
tude has the same weak phase as the b ccs transition. 



S 



sin 2(j)i and A 



are expected if other 



contributions to the decay amplitude can be neglected [J] . 
If, however, the penguin or other contributions are sub- 
stantial, the CP violation parameters for this mode may 
deviate from these values. Employing SU(3) symmetry 
as well as plausible dynamical assumptions, the results 
obtained for B ^ J/tpn^ decay can be used to estimate 
the penguin pollution in B'^ Jhp Kg decay for a very 
precise determination of sin 2(f>i [5| . 

The most recent study of J5° J/^ 7r° decays was 
reported by BaBar [6] using a sample of 232 million BB 
pairs, while the previous Belle analysis [3] was based on 
a data sample corresponding to 152 million BB pairs. 
This measurement of time-dependent CP violation in 
B^ J/^pTT° decays is based on a larger data sam- 
ple that contains 535 million BB pairs, collected with 
the Belle detector at the KEKB asymmetric-energy e+e" 
(3.5 on 8 GeV) colHder Q operating at the T(4S') res- 
onance. The T(45) is produced with a Lorentz boost 
factor of p-f = 0.425 along the 2;-axis, which is anti- 
parallel to the positron beam direction. Since the BB 
pairs are produced nearly at rest in the T(4S') center- 
of-mass system (cms), At is determined from Az, the 
distance between the two B meson decay vertices along 
the z-direction: At = Az/cPj, where c is the speed of 
light. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
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gel threshold Cherenkov counters (ACC) , a barrel-like ar- 
rangement of time-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect 
mesons and to identify muons (KLM). The detector 
is described in detail elsewhere f£|. Two inner detector 
configurations were used. A 2.0 cm radius beam pipe 
and a 3-layer silicon vertex detector were used for the 
first sample of 152 million BB pairs, while a 1.5 cm ra- 
dius beampipe, a 4-layer silicon detector and a small-cell 
inner drift chamber were used to record the remaining 
383 milhon B'B pairs 0. 

We reconstruct J /ij} mesons in the decay chan- 

nel = e or /i) and include up to two bremsstrahlung 
photons that are within 50 mrad of each of the e"*" and 
e~ tracks (denoted as 6+6^(7)). The invariant mass is 
required to be within —0.15 GeV/c^ < Mee(7) — mj/^ < 
-h0.036 GeV/c2 and -0.06 GeV/c^ < M^^ - mj/^ < 
-1-0.036 GeV/c^, where mj/^ denotes the J/4> nominal 
mass [ll| , and Mge(^-) and M^^ are the reconstructed in- 
variant masses from e~^e~{'j) and iJ.'^n~, respectively. 

Photon candidates are selected from clusters of up to 
5x5 crystals in the ECL. Each candidate is required to 
have no associated charged track and a cluster shape that 
is consistent with an electromagnetic shower. To select 
— > 77 decay candidates, the energy of a photon is 
required to be greater than 50 MeV in the ECL barrel and 
100 MeV in the end-cap region. A pair of photons with 
an invariant mass in the range 118 MeV/c^ < M^^ < 150 
MeV/c^ is considered as a 7r° candidate. 

We combine the J/^p and tt*' to form a neutral 
B meson. Signal candidates are identified by two 
kinematic variables defined in the T(4S') rest frame 
(cms): the beam-energy constrained mass Afbc = 
V-^boam ~ iJ^Pi)'^ ^"^^ ^he energy difference AE = 
Y,Ei - Ehcam, where Ehcam = a/s/2 is the cms beam 
energy, and pi and Ei are the cms three momenta and en- 
ergies of the candidate B meson decay products, respec- 
tively. In order to improve the AE resolution, vertex- 
and mass-constrained fits are applied to J/ip 
decays and a mass constrained fit is used for tt*^ 77 
decays. The B meson signal region is defined as 5.27 
GeV/c2 < Mbc < 5.29 GeV/c^ and -0.1 GeV < AE < 
0.05 GeV. The lower bound in AE is chosen to accom- 
modate the negative AE tail of the signal due to shower 
leakage associated with the tt", and to avoid background 
from B^ — > Jli\)K'g (if^ tt^tt") decays. To suppress 
the two-jet-like e^e^ ^ <M (q — u, d, s, c) continuum 
background, we require that the event shape variable, i?2, 
which is the ratio of the second to zeroth Fox- Wolfram 
moment, satisfy R2 < 0.4 [T^ . 

We identify the flavor of the accompanying B me- 
son from inclusive properties of particles that are not 
associated with the reconstructed J/ipn^. The 



algorithm for flavor tagging is described in detail else- 
where |13[. We use two parameters, q defined in Eq.((T]) 
and r, to represent the tagging information. The pa- 
rameter r is an event-by-event Monte Carlo (MC) de- 
termined flavor-tagging quality factor that ranges from 
r = for no flavor discrimination to r = 1 for un- 
ambiguous flavor assignment. It is used only for sort- 
ing data into six intervals. The wrong tag fractions 
for the six r intervals, wi {I — 1,6), and the difference 
in w between and B^ decays. Aw;, are determined 
from data [l3| . The vertex position for the J/^ 7r° de- 
cay is reconstructed using leptons from the J/^ decay. 
The vertex position of /tag is obtained using tracks that 
are not assigned to the J/i/' tt'^ candidate and an inter- 
action point constraint. After all selection criteria are 
applied, we obtain 864 events in the AE — Mbc fit re- 
gion defined as 5.2 GeV/c^ < Mbc < 5.3 GeV/c^ and 
-0.2 GeV < A£; < 0.2 GeV, of which 290 are in the 
signal box. 

We perform an unbinned maximum likelihood fit to the 
AE—Mhc distribution in order to distinguish signal and 
backgrounds. The probability density function (PDF) 
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FIG. 1: Ai5 distribution for events in the A/bc signal region 
(top) and Mbc distribution in the /S.E signal region (bottom) . 
The superimposed curves show the signal (solid line), B 
J/i/j X (dot-dashed line) , combinatorial background (dashed 
line) and the sum of all the contributions (thick solid line). 

of signal is composed of two parts: one is for the can- 
didates that are correctly reconstructed combinations of 
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daughter particles coming from a single neutral B me- 
son, the other corresponds to combinations in which one 
of the final state particles is incorrectly reconstructed 
(i.e. one of the daughter particles originates from the 
other B meson). The former is parameterized by a two- 
dimensional function that is a product of a Crystal Ball 
line shape 



14l | in AE and a Gaussian form in Mbc- This and Anid are fixed to PDG2006 values 



J/i/jKI, other B J/ipX background and combina- 
torial background, respectively. All fractions are func- 
tions of AE and Mbc and are determined from the fit 
discussed above. The PDF for the signal distribution, 
T'sig, is given by Eq. [T] and modified to account for the 
effect of incorrect flavor assignment; the parameters r^o 

Til. The sig- 



parameterization accounts for the fact that AE and Mbc 
distributions are predominantly affected by the shower 
energy leakage in the ECL (in AE) and the beam en- 
ergy spread of the KEKB accelerator (in Mbc). On the 
other hand, the latter is described by a MC-determined 
two-dimensional smooth function. In the signal box, the 
correct combination is estimated to describe 87 ± 2 % of 
the signal events. The background is composed of four 
components: (1) 5° ^ J/^K^, (2) B° J/tpK", (3) 
B J/ipX other than 5° J/iPK°, (4) combina- 
torial background that consists of random combinations 
of particles in BB decays and continuum events. Us- 
ing a large MC sample, the PDFs to describe (1), (2) 
and (3) are determined and then parameterized as two- 
dimensional smooth functions in AE—Al^c- In the fit, we 
fix each yield of the three components, (1), (2) and (3), 
to the values obtained from the MC sample. The dom- 
inant B — > J/ipX contributions, excluding J/i/jKg and 
J/ip K1, come from two-body decays, with well-measured 
branching fractions {B J/ipK^*'>). The combinatorial 
background shapes in AE and Mbc are described by a 
first-order polynomial and an ARGUS function [15i] , re- 
spectively. The purity in the signal region is estimated 
to be 87.9 ± 8.0%. The fractions of J/il^Kl, J/^Kl 
and other J/TpX events are 2.6 ± 0.2%, 2.0 ± 1.2% and 
3.2 ± 0.2%, respectively, while the combinatorial event 
fraction is 4.3 ± 0.5%. The AE and Mbc distributions 
after tagging and vertexing are shown in Fig. [1] 

We determine Sj/^^^o and ^j/^^o by performing an 
unbinned maximum-likelihood fit to the observed At dis- 
tribution: 

JV 

^{Sj/^nO,-^J/4,Tr°) = ]^'^('5j/^TrOi^J/i>7r"; ^^4)7 (3) 
i 

where the product is over all events in the signal region. 
The PDF V is given by, 

V= {l-foi){J d{At')R{AU-At')[ 

fsigPsigiAt ) 

+ fj/iPKs'Pj/rpKs{At') + f.j/^KL'P.I/^KLiAt') 

+fj/^xVj/^xiAt')] 

~\~ ^comb^comb ( ) } 

+fMAU), (4) 

where /sig, fj/^Ks^ Ij/^Kl^ fj/ipx and /comb are the 
fractions of 5° J/V'tt" signal, B° JI^K%, B° 



nal PDF is convolved with the proper-time interval res- 
olution function R{At) 46]. The J/V'^s and 
_B" ^ J/-0 K1 background distributions are described by 
the same Psig, respectively called Vj/^Ks and Vj/^Kl^ 
convolved with R{At). The CP-asymmetry parameters 
Sj/^K%^ -^j/VXg' ^JJAKl and Aj/^^l are fixed to the 
recent Belle results [l7[. The B — > J/ipX background 
excluding the B° J/4>Kl and B^ J/4>Kl com- 
ponents iVj/ipx) is described with an effective lifetime 
as, 
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J/ipX 



{At) 



|At|/r 



J/ii X 



At 



'j/i>x 



{l gAcj,}. 



(5) 



The effective lifetime tj/^x is 1.10 ±0.10 (1.03±0.07)ps 
for the 3(4)-layer-silicon vertex detector sample, which is 
determined by fitting a, B ^ J/ipX MC sample. The 
combinatorial component (Pcomb) is described by a dou- 
ble Gaussian. The relevant parameters are obtained us- 
ing events in the sideband region, 5.20 GeV/c^ < Mbc < 
5.26 GeV/c^ and |A£;| < 0.2 GeV. The fraction /oi and 
PDF Vol describe the outlier component, which is a small 
number of events that have large At values for both sig- 
nal and background. 

The unbinned maximum likelihood fit to the 290 events 
in the signal region results in the CP violation parame- 
ters: 

Sj/i,-nO = -0.65 ± 0.21 (stat) ± 0.05(syst) and 
Aj/^^o = -)-0.08±0.16(stat)±0.05(syst), 

where the systematic uncertainties listed are described 
below. The At distributions and the time-dependent de- 
cay rate raw asymmetry Acp are shown in Fig. [51 where 
Acp = {N+-N_)/{N+ + N_) and N+ (iV_) is the num- 
ber of candidate events with q = +1 (—1). 

The systematic errors are fisted in Table H The mam 
contributions to the systematic error in Sj/^^^o are due 
to uncertainties in the vertex reconstruction and to a 
small fit bias. The vertex reconstruction systematic er- 
ror consists of uncertainties in the interaction point pro- 
file, charged track selection based on the track helix er- 
ror, helix parameter corrections, event selection based on 
At and goodness of fit in the vertex reconstruction, and 
the small SVD misalignment. The systematic uncertain- 
ties due to the parameters wi and Awi are estimated by 
varying the parameters by their one standard deviation 
(cr) errors. We vary each resolution function parameter 
by ±1(7 and assign a systematic error as the quadratic 
sum of the resulting deviations in S and A. The fit bias 
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5 At(ps) 



FIG. 2: At distribution of B° 
q = +1 (a) and q — —1 (b). 



~» J/ip 7r° candidate events for 
The dashed Unes are the sum 
of backgrounds while the sohd hnes are the sum of signal and 
backgrounds, (c) is the raw asymmetry (Acp) distribution. 
The curve is the projection of the fit result. 



TABLE I: Systematic uncertainties 



Parameter 






Vertexing 


±0.050 


±0.034 


Wrong tag fraction 


±0.009 


±0.009 


Resolution function 


±0.008 


±0.007 


Fit bias 


±0.013 


±0.010 


Physics parameters 


±0.004 


±0.001 


B ^ J/ip X CP asymmetry 


±0.004 


±0.001 


PDF Shape and fraction 


±0.009 


±0.005 


Background At shape 


±0.006 


±0.001 


Tag side interference 


±0.001 


±0.038 


Total 


±0.054 


±0.054 



systematic error is evaluated from an ensemble of MC 
samples as the difference between the input and fitted 
values of S and A. The errors in the physics parameters 
Tgo and Arud are taken into account. To estimate the 
systematics from i?° J/ip Kg and ^ J/tp , we 
vary their fractions and CP asymmetry parameters, S 
and A by ±lcr. We estimate the systematic uncertainty 
from the B J/tp X backgrounds other than J/ -0 by 
scaling the systematic errors due to J/ -0 and J/V' 
according to the amount of background contamination. 
The Ai? and Afbc parameters and fractions of signal and 
backgrounds are varied to estimate the systematic errors. 
We vary the MC-determined parameters by ±2cr to take 
account of possible imperfect modeling in MC. We in- 
clude the effect of tag side interference ;1S]. The tag 



side interference is caused by the interference between 
the two amplitudes of B decays into charmed mesons, 
i.e. caused by Vch and V^h- Therefore it is expressed 
by four parameters. Tint (size of interference between Vch 
and Viib amplitudes), 0i, 03 and 5 (strong phase differ- 
ence between Vcb and V-ah mediated amplitudes). Since 
this interference results in a potential direct CP viola- 
tion, Aj/^Tr° is much more affected than Sj/^t^o. We 
sum each of the contributions in quadrature to obtain 
the total systematic error. 




-1 -0.8 -0.6 -0.4 -0.2 
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FIG. 3: The confidence regions for Sj/^^o and Aj/^^o. The 
contours correspond to 1-C.L. = 3.17 x 10~^ (lo")) 4.55 x 
10"^ {2a), and 2.70 x 10"^ {3a). The point with error bars 



corresponds to the measured S 



J /in 



and A 



circle is a boundary derived from Eg.©. 



j/ii 1 



values. The 



The confidence regions of our measurement in the 
Sj/iPttO and A,j /^t^o plane are shown in Fig. [3] We 
evaluate the statistical significance of this CP asymme- 
try measurement using a two-dimensional Fcldman and 
Cousins method [l3|, taking both statistical and sys- 
tematic uncertainties into account. We found that our 
Sj/jp^o measurement has a significance greater than 2.4 
tJ for any Aj /^t^o value. 

In summary, we measure the CP violation parame- 
ters in — > J/-07r° decays using 535 x Kf'BB pairs: 



S 



-0.65 ± 0.21 (stat) ± 0.05(syst) and Aj/^^a 



±0.08 ± 0.16(stat) ± 0.05(syst). We measure mixing- 
induced CP violation with 2.4 a significance. This re- 
sult supersedes our previous measurement Q and ex- 
hibits significant improvement in precision compared to 
the latest BaBar measurement It is consistent with 
the measured value of sin 20i in 6 ^ ccs decays 17 , lO] , 
as expected in the Standard Model. 
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